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Abstract 

The factors regulating the structure of food webs are a central focus of community and 
ecosystem ecology, as trophic interactions among species have important impacts on 
nutrient storage and cycling in many ecosystems. For soil invertebrates in grassland 
ecosystems in the Netherlands, the site-specific slopes of the faunal biomass to organism 
body mass relationships reflected basic biochemical and biogeochemical processes 
associated with soil acidity and soil C:N:P stoichiometry. That is, the higher the 
phosphorus availability in the soil, the higher, on average, the slope of the faunal 
biomass size spectrum (i.e., the higher the biomass of large-bodied invertebrates relative 
to the biomass of small invertebrates). While other factors may also be involved, these 
results are consistent with the growth rate hypothesis from biological stoichiometry that 
relates phosphorus demands to ribosomal RNA and protein production. Thus our data 
represent the first time that ecosystem phosphorus availability has been associated with 
allometry in soil food webs (supporting information available online). Our results have 
broad implications, as soil invertebrates of different size have different effects on soil 
processes. 

Keywords: allometry, biological stoichiometry, biomass-size spectra, carbon, land use change, mass- 
abundance relationships, nitrogen, phosphorus, soil food webs 

Received 15 December 2008; revised version received 14 February 2009 and accepted 16 February 2009 



Introduction 

Biological stoichiometry (Elser et ah, 1996, 2000a, b; Ster- 
ner & Elser, 2002) is the study of the balance of energy 
and multiple chemical elements in living systems. A 
central idea of biological stoichiometry is the growth 
rate hypothesis (GRH), which proposes that variation in 
C : N : P ratios among organisms reflects differential allo- 
cation to P-rich ribosomal RNA in support of rapid 
growth rate. Thus, fast-growing taxa have low C:P 
and N : P biomass ratios, making them more susceptible 
to P-based stoichiometric food quality constraints (Elser 
et al, 2000b; Sterner & Elser, 2002). As a result, ecosystem 
conditions that produce organic matter with high C:P 
and N : P ratios are likely to result in inefficient trophic 
transfer, reduced biomass of upper trophic levels, and 
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potential deterministic extinction of P-rich herbivores 
(Perez-Moreno & Read, 2001; Sterner & Elser, 2002). 
Many of these ideas have been extensively explored in 
aquatic systems (Blanco et al., 1998; Sterner & Schulz, 
1998; Gamble et al, 2006; Eyto & Irvine, 2007). 

Body size relations are important for understanding 
ecological processes (Peters, 1983; Calder, 1984; Mulder, 
2006; Damuth, 2007), but ecological effects of C : P and 
N:P stoichiometric imbalance have not been widely 
considered for the differently sized soil invertebrates 
that make up belowground food webs. In aquatic 
habitats, in fact, allometric associations of nutrient con- 
tent with abundance and biomass have already been 
established (e.g., Sheldon et al, 1972; Blanco et al, 1998). 
In contrast, in terrestrial ecosystems the connections 
between nutrient stoichiometry and size scaling are 
only now beginning to be considered (Schade et al, 
2003; Kerkhoff et al, 2005; Meehan, 2006; Enquist et al, 
2007; Martinson et al, 2008; Reuman et al, 2009). 
Eurthermore, the size spectra of different clades in 
terrestrial ecosystems have traditionally been studied 
separately, often from an entomological perspective 
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(Siemann et ah, 1999), obscuring community- wide allo- 
metric patterns. Here we investigate if the distribution 
of biomass across body-size classes of all occurring soil 
invertebrates (henceforth, the faunal biomass spectrum) 
and their microbial resources are associated with nu- 
trient concentrations and soil C : N : P ratios. In P-defi- 
cient sites, bacteria are expected to grow quite slowly 
and thus to have reduced P content due to lower RNA 
concentrations (Elser et ah, 2000a, 2003; Gillooly et al., 
2005), potentially imposing stoichiometric food quality 
constraints and impairing the development of higher 
trophic levels. The numerical abundance of bacteria 
reflects abiotic conditions and is positively correlated 
with soil pH (Mulder et al., 2005a, b). To evaluate the 
possible impacts of soil conditions on the structure of 
soil food webs we examined the relationships between 
soil pH and C:N:P ratios (in mass units) and the 
regression slopes of both the faunal biomass spectra 
(which plot log biomass by body-mass categories as a 
function of log body mass) and the mass-abundance 
planes (which plot log numerical abundance as a func- 
tion of log body mass) for a large number of grassland 
soils encompassing broad ecological gradients in the 
Netherlands. 

Materials and methods 

Sites and organisms 

One recent study showed that grasslands under organic 
management might exhibit a departure from power law 
behaviour due to organic fertilizer inputs (Reuman 
et al., 2008). To avoid such effects, we selected, from 
all the locations monitored by the RIVM (Bilthoven, the 
Netherlands), 12 of the organic grasslands under low- 
intensive, bio-dynamic management and 10 ex-organic 
farms that were abandoned for at least a decade (Table 
SI). In contrast to soil acidity and P, which typically 
decrease following land abandonment (Fig. 1), soil 
nutrient ratios are expected to increase after abandon- 
ment (Knops & Tilman, 2000). Environmental gradients 
were chosen to obtain a certain overlap in all the 
investigated abiotic correlates (Table 1). 

Eor this study, upper horizons (top 10 cm) of the 22 
Pleistocene sandy soils selected across the Netherlands 
were randomly sampled between 2002 and 2005. Hor- 
izons have similar bulk density (1 .20pb i 0-29 SD) and 
organic matter content (7.78% ± 3.69 SD). Soil acidity of 
oven-dried soil samples was measured in 1 M potassium 
chloride solution (Table SI). Soils were analysed, at least 
in triplicate, after preliminary wet digestion. Ntot and 
Ptot (henceforth N and P) were determined by colori- 
metric analysis. Soil P content was also measured as 
phosphate (Pw) by using aqueous extraction at water-to- 



CD 



(J) 




Abandoned grassland 



Fig. 1 The development of mature grasslands is a primary con- 
cern of many restoration efforts in the Netherlands. We investi- 
gated both managed grasslands (average soil pH 5.3 and C : N : P 
ratio 41.1 :2.8: 1 in mass units) as abandoned grasslands (average 
pH 4.5 and C : N : P ratio 98 : 5.2 : 1 in mass units). Holcus lanatus, 
Lolium perenne, and Dadylis glomemta characterise the studied 
grasslands together with Alopecurus geniculatus and Poa trivialis. 

Table 1 Soil C : N : P stoichiometry and abiotic properties 



Main abiotic properties Abandoned Managed 



Carbon (gkgsoildry mass~^) 


54.0 




23.7 


37.9 


± 


17.1 


Nitrogen (gkgsoildry mass~^) 


2.9 




1.0 


2.5 


± 


0.7 


Phosphorus (g kg soil dry mass" ^) 


0.6 




0.3 


0.9 


± 


0.4 


Carbon-to-phosphorus ratio 


98.0 




58.4 


41.1 




8.8 


Carbon-to-nitrogen ratio 


18.5 




4.7 


14.7 




2.3 


Nitrogen-to-phosphorus ratio 


5.2 




2.2 


2.8 




0.7 


Soil organic matter (%) 


9.3 




4.1 


6.5 




2.9 


Soil acidity (pH in KCl) 


4.5 




0.3 


5.3 




0.4 



Means ± SE for 10 abandoned grasslands (mature meadows) 
and 12 biomanaged grasslands under organic regime (com- 
plete data online in Table SI). 



soil ratio of 60 : 1 by volume, after 22 h of pre-equilibra- 
tion with soil with water and 1 h of gentle shaking 
before filtration (Sissingh, 1971). Average values of C 
(dry combustion) were obtained by conversion of soil 
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Table 2 Faunal web structure, soil microbiology and biodi- 
versity distribution 



Main biotic properties 



Abandoned Managed 



Total number of potential trophic 

links (linkages between genera) 
Biodiversity of nematodes 

(Nematoda) 
Biodiversity of soil mites (Acarina) 
Biodiversity of springtails (Insecta) 
Biodiversity of enchytraeids 

(Oligochaeta) 
Fungal-to-bacterial ratio 

(^gC^ig-^C) 



710 ± 159 706 ± 225 



28 ±3 

16 ±3 
6±3 
6 ± 1 



25 ±5 

16 ±5 
10 ±3 
5±1 



0.15 ± 0.08 0.26 ± 0.31 



Means ± SE for 10 abandoned grasslands (mature meadows) 
and 12 biomanaged grasslands under organic regime (com- 
plete data online in Table S2). 

organic matter (SOM) using SOM = Corg x 1.724 (the 
traditional Van Bemmelen conversion factor). 

The total abundance of bacterial cells was determined 
in duplicate by fluorescent staining [5-(4, 6-dichloro- 
triazin-2-yl) aminofluorescein] and combined micro- 
scopy and automatic image analysis. For the 
conversion of bacterial cell volume (}im^) to dry bio- 
mass (|ig) v^e used the biovolume-to-carbon factor of 
Van Veen & Paul (1979), assuming a bacterial carbon 
content of 50% (Herbert, 1976; Mulder et al, 2005a; Table 
2). Hyphae v^ere viev^ed with fluorescent staining and 
direct microscopy at x 2500 and the total hyphal length 
was estimated by the line intercept method. The myce- 
lium was inferred from the length of the counted 
branches, assuming a hyphal diameter of 2.5 \mv (Table 
S2). All fungal measurements used direct microscopy 
and palynological treatments (Mulder et al., 2005b). 

The soil invertebrates (Table S2) were measured as 
follows. For worm-like animals (nematodes and enchy- 
traeids), the body length and width of at least 150 
nematodes were measured to the nearest 5 |im with an 
eyepiece micrometer; all enchytraeids were measured 
individually, including juveniles and resting stages. The 
soil nematode samples were collected randomly across 
the investigated site. In each of the 22 locations, a 
composite sample was obtained by mixing 320 soil cores 
(diameter 2.3 cm x 10 cm) in a plastic container, and 
approximately 500 g of wet soil was collected in glass 
jars. Nematodes were extracted, within 1 week, from 
100 g of wet soil using the funnel elutriation comple- 
mented with sieving and cottonwood extraction after 
crawling through a cotton filter over a 2-day period. Two 
clean suspensions in 10 mL water were screened with a 
stereoscope to count the individuals; the total of nema- 
todes was estimated by counting 10% of the extracted 
animals twice. Permanent mounts in 4% formaldehyde 



were made on mass slides. Nematodes were identified 
by light microscopy at x 400-600 (Mulder et al, 2003, 
2005b). Microarthropods were extracted from the soil by 
placing six discs of the soil sample in a TuUgren funnel 
(Siepel & Van de Bund, 1988; Rombke et al, 2006). The 
temperature in the upper part of the funnel was set at 
30 °C and kept at 5 °C in the lower part. The organisms 
moved downwards to escape the heat, dropped through 
a funnel and were collected into a bottle containing 70% 
ethanol. The total extraction time was 1 week. For each 
sample, 70 individuals were counted and identified at 
X 200-1000 via a gel-based subsampling method (Jagers 
op Akkerhuis et al, 2008). Extracted arthropods were 
divided in body-size classes to estimate the correspond- 
ing individual body mass (M in microgram of dry 
weight). All values were averaged over the life-stages 
of all the counted individuals AT (Mulder et al, 2005a, b). 
These body-mass averages at genus level offer the best 
available combination of high environmental informa- 
tion and low noise (Mulder et al, 2006, 2008). 

Allometric scaling 

Descriptors derived from allometric patterns require a 
brief explanation. If the species or genera of a local 
community are plotted as points in a plane with abscis- 
sa (horizontal axis) log(M), and with ordinate log(N), 
then the points have been found to fall approximately 
along a straight line with negative slope. This mass- 
abundance slope is the coefficient a in the linear model 
log(N) =ax log(M) + h fitted to data from a single soil 
community food web. The data can comprise organisms 
belonging to different kingdoms or other taxonomic 
groups. We computed the mass-abundance slopes for 
all soil invertebrates (A, from Animalia), and added two 
further aggregated points, one for Fungi (F, resulting in 
a scatter of all A + F points) and one for Eubacteria (E, 
resulting in a A + E + F scatter of all recovered soil 
eukaryotes and bacterial cells). Because the biomass of 
a taxon is its numerical abundance times its body mass, 
B = NM, and because increasing M is associated with 
increasing trophic height (Mulder et al, 2005a, 2008; 
Reuman et al, 2008), the mass-abundance slopes for A, 
A + F, and A + E + F are supposed to indicate how 
(faunal) biomass and population density change for 
(faunal) taxa of increasing trophic height. Mass-abun- 
dance relationships reveal how the size structure of a 
community food web interacts with its trophic structure 
and biomass distribution (White et al, 2007). If the 
resulting mass-abundance slope is exactly —1, then 
the trend is for all taxa to have equal biomass. All 
identified genera (Table 2) fell into 22 trophic guilds. 
Trophic links among guilds (from resource to consu- 
mer) were inferred from published literature (Table S2). 
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It was assumed that every taxon in a resource guild was 
trophically linked to every taxon in a consumer guild. 

Merging the classic allometric formula log(N) 
= ax log(M) + b with log(B) = log(M) + log(N), we ob- 
tain log(B) = log(M) + ax log(M) + b = {l+a)x log(M) 
+ b (e.g., Cyr et al, 1997; Mulder et al, 2008) Theoreti- 
cally, thus, the slope of a faunal biomass spectrum is 1 
plus the slope of the relation between log(N) and log(M) 
for all the occurring animals; e.g., a zero slope of the 
biomass spectrum (constant biomass across trophic 
levels) corresponds to a slope of —1 in a plot of log(N) 
in bins as a function of log(M) (Rossberg et al., 2008). To 
compute the faunal biomass spectrum slopes using bins 
of constant linear width, the range from the smallest 
log(M) to the largest log(M) occurring in any of the 
grasslands was divided into 10 equal bins. For each bin, 
the log of the total faunal biomass was computed. All 
values for bins containing fauna (bins with zero obser- 
vations are excluded) were regressed against the log(M) 
at the centre of each bin. To test whether binning 
introduces a bias towards small or large organisms as 
suggested by White et al. (2008), we compared (binned) 
biomass spectrum slopes with three (unbinned) mass- 
abundance slopes. Computations used SAS version 9.1.3, 
PC-ORD version 4.20 and the excel Visual Basic optimiza- 
tion toolbox. We set 1 % significance to detect violations. 



P becomes curtailed under lower soil pH (Bohn et al., 
1985), as in our abandoned grasslands. 

Soil acidity of abandoned grasslands was almost 
seven times that of managed grasslands (pH averaged 
4.48 ± 0.28 and 5.32 ± 0.36 SD, respectively). Soil pH 
contributes strongly to the explanation of the allometric 
variation within a pH range of 3.8-6.0 (Fig. 2 and Table 
SI). We tested the resulting correlation between soil pH 
and the biomass spectrum slope by adding 10 locations 
belonging to a completely different ecosystem, namely 
dry heathlands with very low soil pH (3.15 ± 0.23 SD), 
to extend the investigated pH range down to 2.9. We 
obtained for this recalculated linear regression a slope 
of 0.375 ± 0.025 SF, which is statistically indistinguish- 
able from the original slope of 0.351 ± 0.046 SF with a 
99% CI. 

Managed grasslands supported invertebrate commu- 
nities in which biomass increased with the central 
log(M) values of bins in biomass spectra. Steep faunal 
biomass slopes were associated here with high popula- 
tion densities of enchytraeids that fall at the upper end 
of the range of body sizes (Table S2). In contrast, 
abandoned grasslands supported invertebrate commu- 
nities in which biomass increased less with increasing 
central log(M) values of bins in biomass spectra. Shal- 
low faunal biomass slopes were associated with higher 



Results and discussion 

There were significant differences in the concentrations 
of soil macronutrients. The C and N concentrations 
were closely correlated (Table SI) and highest, on 
average, in the abandoned grasslands (i.e., not fertilized 
for at least one decade) and lowest, on average, in 
the managed grasslands (Table 1). Total soil P contents, 
on the contrary, were highest in the managed grass- 
lands, due to inputs of organic fertilizers, and lowest 
in the abandoned grasslands (939 ± 400 vs. 627 ± 
277SDmgPkg"^). Coefficients of variation (100 x SD/ 
mean) of total soil P were 43% and 44% for managed 
and abandoned grasslands, respectively. Soil P content 
was uncorrelated with soil pH, in contrast to the phos- 
phate content after water extraction [Pearson's correla- 
tion coefficients 0.024 for P and pH (P>0.05) and 0.53 
for Pw and pH (P = 0.0115), respectively]. Phosphate 
content after water extraction seemed to reflect the 
maximal possible concentration of phosphorus in the 
thin layer of water around soil particles (35.3 ± 9.1 vs. 
19.6 ± 24.6 SDmgL"^ in managed and abandoned 
grasslands, respectively). Biologically available P is 
thought to increase within a pH range of 5-6 (Chapin 
& Fviner, 2003; Sims & Pierzynski, 2005; Cleveland & 
Liptzin, 2007), as in our managed grasslands where 
fertilization and liming occurred regularly. In contrast. 
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Fig. 2 The slope of the biomass spectrum of soil invertebrates 
as predicted by soil acidity (pH in KCl). We combined the body 
masses of all individuals whose size fell within an equal interval 
of body size and plotted, separately for each grassland, the log 
summed biomass of all animal taxa within log(M) bins against 
bin centres. In this way, the slopes of the linear regression of log 
summed biomass become a function of the bin centres on a 
log(M) scale. Pearson's correlation coefficients show a significant 
correlation only between the soil pH and the logarithm of C : P 
(P = 0.0010). Managed grasslands are squares and abandoned 
grasslands are diamonds. 
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Fig. 3 Different allometric scalings and soil nutrient ratios. Arranged according to a decreasing coefficient of variation in our sites 
(slopes rescaled in degrees, not shown), from top to bottom: faunal biomass spectrum slope (A = Animalia), faunal mass - abundance 
slope (A = Animalia), eukaryotes' mass - abundance slope (A + F = Animalia and Fungi), and mass - abundance slope of the complete 
community food web (A + E + F = Animalia, Eubacteria, and Fungi). The logarithm of the C : P ratio equals the sum of the log C : N and 
log N : P ratios and is the best sole predictor for allometric scaling after pH. Symbols as in Fig. 2. 



population densities of nematodes that fall at the lower 
end of the size range, indicating that biomass within 
log(M) bins increased with average body size but 
increased less rapidly on average than after application 
of manure and lime. Similar associations of soil pH with 
different taxa and communities have been recorded 
recently (Mulder et al, 2003, 2005b; Fierer & Jackson, 
2006). 



Plotting biomass spectrum slopes and mass-abun- 
dance slopes against the logarithms of the soil nutrient 
ratios revealed strong relationships between nutrient 
ratios and faunal biomass size spectra (Fig. 3). The 
interpretation of the biomass spectrum slope is intui- 
tive: negative trends indicate that faunal biomass under 
nutrient deficiency declined with increased body size. 
Biomass spectrum slopes were most significantly corre- 
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lated with the logarithm of the soil C : P ratio (Pearson's 
correlation coefficient -0.7013, P = 0.0003), but less with 
the logarithms of the N : P and C:N ratios (P = 0.0049 
and P = 0.0099, respectively). As a matter of fact, in all 
cases of allometric scaling (n = 22, four kinds of lump- 
ing, resulting in 88 investigated cases comprehending 
two ecosystem types), the logarithm of the C:P ratio 
was the best sole predictor (Fig. 3, right column). 

For further interpretation of mass-abundance slopes, 
derived from unlumped and unbinned points with 
specific numerical abundance and body-mass average, 
one brief example may be useful. Let us imagine a 
simple food web consisting of four genera, namely 
one nematode, one mite, one springtail, and one enchy- 
traeid. Let their respective dry weights be 0.1, 1, 10, and 
100 )ig on average (extensive empirical data in Table S2). 
After log- transformation, their log(M) will become — 1, 
0, 1, and 2. Seen that in a plane with abscissa log(M) and 
with ordinate log(N) most populations have been found 
to fall approximately along a straight line with negative 
slope, we assume for ease of computation that the popu- 
lation densities of these four genera are equal to 10000, 
1000, 100, and 10 individuals per square meter, respec- 
tively. After log-transformation, their log(N) will become 
4, 3, 2, and 1. Their specific log(B) equals log(M) + 
log(N)=-l +4 = 0 + 3 = 1 +2 = 2 + 1 =3. Hence, these 
four populations wiU keep a dry biomass of 1 mg m~^ 
and, if plotted in a log(M)— log(N) plane, the result- 
ing linear regression slope will be exactly equal to — 1. 
Our mass-abundance slopes were always negative 
(— 1<^7<0), as expected from the positive biomass spec- 
trum slope. The question which arises here is: how do 
slopes derived from such mass-abundance relationships 
reflect soil resource limitation? 

The linear regressions of log(N) as function of log(M) 
account for nearly 50% of the variation in log(N) of 
organisms in most grasslands, but linear regressions 
may account for only 30% of the variation in some 
managed grasslands. The expected numerical abun- 
dances of ubiquitous smallest animals (nematodes) 
vary more than the numerical abundances of the larger 
animals (mites, springtails, and enchytraeids; see Table 
S2). Numerical abundances of smaller soil invertebrates 
differ up to four orders of magnitude from those of 
larger invertebrates. Moreover, the body-mass ranges 
of soil invertebrates reflect dominant feeding habits, 
since opposite functional types coming out of taxo- 
nomic inventories can be recognized: 55% of the 
microarthropods graze, browse or ingest fungal re- 
mains and 40% of nonparasitic nematodes graze on 
bacterial cells (Mulder et al, 2005b and Table S2). 

Therefore the allometric scaling strongly suggests 
that the microfauna (nematode populations) copes 
better with P-limited conditions than the mesofauna 



(microarthropod populations). Equivalently, the less the 
limitation by either nitrogen or carbon relative to phos- 
phorus, or the greater the phosphorus limitation, the 
more rapidly the population density of invertebrate 
genera decreases as the average body mass of the 
invertebrate genera increases. All the negative trends 
along the environmental correlates of Fig. 3 indicate 
thus (i) that faunal biomass declined with increased 
body size within log(M) bins, and (ii) that faunal 
population density decreased with increased body size. 
Hence differences in soil P availability are reflected 
within our grasslands (and possibly in other biomes 
as well) in the body-size distribution of the biomass and 
density of the soil microfauna and mesofauna, since 
stoichiometric theory predicts that soil fauna with high- 
er P demands would suffer a competitive disadvantage 
in lower P soils due to poorer stoichiometric food 
quality The GRH (Elser et al, 2000a, 2003) suggests that 
high P content and low C:P ratio in biomass reflect 
increased allocation to P-rich ribosomal RNA, which in 
turn enables increased protein synthesis, increased 
growth rates of individuals, and increased population 
growth rates. Empirical investigations support the GRH 
in varied contexts (e.g., Elser et ah, 2000a, 2003). The 
GRH provides a natural context for interpretation of at 
least some of our observations. Agricultural soils with 
higher P content, as a result of edaphic conditions, land 
use and history (in our case, increased P availability due 
to fertilization and liming effects on pH), have a twofold 
lower average soil C:P ratio than abandoned grass- 
lands (Table 1). 

As recently stated by Urabe & Waki (2009, p. 529), 
'herbivores tend to increase the ingestion rate when 
food with low nutrient content relative to C is supplied'. 
Many larger-bodied soil herbivores are further known 
for a certain combination of intrinsic and microbial 
sources of cellulolysis in their diet (e.g., Douglas, 2009, 
p. 42), and most growth and reproduction rates have 
been observed to decline as a function of body mass 
(Hendriks & Mulder, 2008; Beardall et al, 2009). We 
suggest that increase in P availability relaxes stoichio- 
metric food quality constraints on small-bodied pri- 
mary consumers (here, microbivores), which increases 
trophic transfer efficiency, and supports larger-bodied 
consumers (including herbivores) also at higher trophic 
levels. That is, soils richer in the key nutrient P have 
organic matter that is of higher quality for invertebrates, 
allowing the maintenance of food webs with increased 
abundances of larger taxa. As a result, biomass spectra 
in P-rich soils have steep positive slopes while relatively 
inefficient food webs in stoichiometrically imbalanced 
low-P soils support mostly small-bodied fauna and thus 
produce shallow (and theoretically even negative) bio- 
mass spectrum slopes. The extent to which low-P soils 
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support the population density of the small-bodied 
fauna is shown by the more negative mass-abundance 
slopes. 

Conclusions 

Alterations in soil C : N : P ratios, such as those driven 
by direct anthropogenic influences like nitrogen deposi- 
tion, liming, and fertilizer application (Hunt & Wall, 
2002; Wardle, 2002; Mulder et al, 2005b, 2008; Persson 
et al., 2008) and /or those driven by climatic changes like 
increased atmospheric CO2, global warming and tor- 
rential rainfall (Wardle et al, 1998; Hunt & Wall, 2002; 
Voigt et al, 2003, 2007; Sardans & Penuelas, 2007; Urabe 
& Waki, 2009), may affect the rates at which soil biota 
carry out ecosystem services by affecting the faunal 
biomass distribution. However, these ideas require 
further investigation and we require more data on the 
C:N:P stoichiometry, threshold elemental ratios, 
growth rates, and RNA demands of soil fauna. 

The association of the faunal biomass spectrum to soil 
acidity is more than a chemical reaction to a relative 
concentration of [H^]-ions since soil pH has a strong 
impact on the nutrient mobility, adsorption, and pre- 
cipitation. Our 10 abandoned grasslands (lower soil pH) 
have a higher concentration of [H^]-ions than grass- 
lands managed with lime, in which the cation exchange 
capacity increases. Furthermore, our data are correla- 
tive in nature. Thus, the hypothesis that stoichiometric 
imbalance imposed by P-limitation alters soil biomass 
spectra is in need of experimental test. Such a test 
would involve quantifying biomass spectra in experi- 
mental plots subjected to different levels of nutrient 
enrichment. Such experiments seem particularly timely 
given the recent indication that P-limitation of primary 
production may be more widespread in terrestrial eco- 
systems than previously appreciated (Elser et al, 2007). 

In summary, faunal biomass spectrum slopes were 
strongly associated with liming practice, soil nutrient 
conditions (as reflected by phosphate), and organic 
fertilizers (as reflected by the soil C : P ratio). Previous 
studies on faunal biomass size spectra in aquatic eco- 
systems reported lower slopes in systems subjected to 
heavy fishing and hypothesized that fishing decreased 
the slope of a biomass spectrum by preferentially re- 
moving larger organisms (Pope et al, 1988; Jennings, 
2005), a 'top down' explanation. Our data showing that 
soil P availability is associated with increasingly posi- 
tive slopes of the faunal biomass size spectrum in soil 
systems represent the first observation from soil food 
webs for a l3ottom up' effect on soil food-web size 
structure. The GRH, with or without increased trophic 
transfer efficiencies, provides a natural interpretation of 
why (artificially) P-enriched systems have food webs in 



which biomass increases with consumer body size 
while biomass declines with body size in P-limited 
habitats. These insights potentially connect biogeo- 
chemistry with the allometric relationships between 
biomass, abundance, and body size. 
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